Prediction and control of elevated temperatures within landfills under aeration and recirculation based on the thermal non-equilibrium model.

Aeration is an effective approach to sustainable landfilling but may lead to elevated temperatures within landfills, resulting in landfill fires or explosions. Therefore, aeration is usually combined with leachate recirculation to control the elevated temperatures within landfills. To predict landfill temperatures during aeration and recirculation, a local thermal non-equilibrium model is developed considering the heat generation of biodegradation, the heat removal due to evaporation and leachate-gas flow, and the effects of the capillary. The solver is implemented in OpenFOAM based on the finite volume method and validated against a waste-column experiment and an in-situ aeration test. The simulation results demonstrate that the assumption of local thermal equilibrium will distinctly overestimate the temperature, maximally by 15 °C in the studied case. The model is then used to simulate a typical aerobic landfill unit to investigate the formation of explosive gas mixtures and elevated temperatures under different operating conditions. The simulation results of gas composition suggest that aeration will not result in explosive gas within landfills. A reasonable recirculation method for temperature control with corresponding operating parameters under a group of values of aeration pressure (2000-4000 Pa) and recirculation rate (0.0001-0.0008 m/s) are proposed, which can provide some guides for the design of an aeration and recirculation combined system. For a given total volume of added leachate, a higher recirculation rate does not always mean better cooling, and the cooling effect of continuous recirculation is better than that of intermittent recirculation.

Modeling the spreading and remediation efficiency of slow-release oxidants in a fractured and contaminated low-permeability stratum.

Traditional remediation technologies cannot well remediate the low permeability contaminated stratums due to the limitation in the transport capacity of solute. The technology that integrates the fracturing and/or slow-released oxidants can be a new alternative, and its remediation efficiency remains unknown. In this study, an explicit dissolution-diffusion solution for the oxidants in control release beads (CRBs) was developed to describe the time-varying release of oxidants. Together with advection, diffusion, dispersion and the reactions with oxidants and natural oxidants, a two-dimensional axisymmetric model of solute transport in a fracture-soil matrix system was established to compare the removal efficiencies of CRB oxidants and liquid oxidants and to identify the main factors that can significantly affect the remediation of fractured low-permeability matrix. The results show that CRB oxidants can achieve a more effective remediation than liquid oxidants under the same condition due to the more uniform distribution of oxidants in the fracture and hence a higher utilization rate. Increasing the dose of the embedded oxidants can benefit the remediation to some extent, while at small doses the release time over 20 d has little impact. For extremely low-permeability contaminated stratums, the remediation effect can be significantly improved if the average permeability of the fractured soil can be enhanced to more than 10-7 m/s. Increasing the injection pressure at a single fracture during the treatment can enlarge the influence distance of the slow-released oxidants above the fracture (e.g., 0.3-0.9 m in this study) rather than below the fracture (e.g., 0.3 m in this study). In general, this work is expected to provide some meaningful guidance for the design of fracturing and remediating low permeability contaminated stratums.

Influence mechanism of thermally enhanced phase change on heat transfer and soil vapour extraction.

Thermal enhanced soil vapour extraction (T-SVE) is a remedial technique involving gas, aqueous, solid and nonaqueous phases along with mass and heat transfer. Interphase mass transfer of contaminants and water evaporation/condensation will cause the redistribution of phase saturation, eventually affecting the performance of T-SVE. In this study, a multiphase, multicompositional and nonisothermal model was developed to simulate the T-SVE treatment of contaminated soil. The model was calibrated using published data from the SVE laboratory and T-SVE field experiments. The temporal and spatial distributions of the contaminant concentrations in the four different phases, the mass transfer rates and the temperatures are presented to reveal the coupling interactions that occur between multiple fields during T-SVE. A series of parametric studies were carried out to investigate the effect of water evaporation and adsorbed/dissolved contaminants on the T-SVE performance. It was found that endothermic evaporation, exothermic condensation and the interaction between different removal paths of a contaminant played critical roles in the thermal enhancement of SVE. Ignoring them can result in significant differences in the removal efficiency values.

Forward and back diffusion of reactive contaminants through multi-layer low permeability sediments.

Contaminants stored in the low permeability sediments will continue to threaten the adjacent shallow groundwater system after the aquifer is remediated. Understanding the storage and discharge behavior of contaminants in the aquitards is essential for the efficient remediation of contaminated sites, but most of the previous analytical studies focused on nonreactive solutes in a single homogenous aquitard. This study presents novel analytical solutions for the forward and back diffusion of contaminants through multi-layer low permeability sediments considering abiotic and biotic environmental degradation. Three representative source depletion patterns (i.e., instantaneous, linear, and exponential patterns) were selected to describe the dissolution of dense non-aqueous phase liquids (DNAPL) in the aquifer more realistically. At the forward diffusion stage, the mass storage of contaminants in the aquitards with the instantaneous pattern is the largest, nearly twice that with the exponential pattern. A simple equivalent homogeneous model is generally adopted in the risk assessment. However, relative to the proposed multi-layer model, it will significantly underestimate the onset of the back-diffusion of heterogeneous aquitards and overestimate the persistence of aquifer plumes. The previously-reported semi-infinite boundary assumption is also not applicable, with a maximum error of over 200% in the long-term prediction of back diffusion behavior of a thin aquitard. Moreover, when the degradation half-life is less than 16 years, less than 10% of the contaminants stored in the aquitards will diffuse into the overlying aquifer, suggesting that biostimulation or bioaugmentation can effectively mitigate back-diffusion risk. Overall, the proposed diffusion-reaction coupled model with multi-layer media is of great value and high demand in predicting the back-diffusion behavior of heterogeneous aquitards and guiding the soil bioremediation.

Enhanced delivery of amendments in fractured clay sites based on multi-point injection: An analytical study.

Fracturing technology that can enhance the delivery of amendments has attracted attention in the remediation of low-permeability contaminated sites. However, there are few works on the enhanced delivery of amendments based on multi-point injection in a fracture-matrix system. This study develops a two-dimensional analytical model for enhanced delivery of amendments in a finite-domain low-permeability matrix through multi-point injection in a natural, hydraulic or pneumatic fracture. The mechanisms of advection, diffusion, dispersion, sorption and degradation are considered in the model and any injection form (e.g., pulse injection, periodic injection or slow-release injection) can be embedded to obtain a specified solution. Then, a new linear factor R*, which is the ratio of the peak concentration to the trough concentration on the same plane, is introduced to evaluate the concentration fluctuation in the fracture and matrix. Results show that with a stronger line source formed in the fracture right after injection (corresponding to a small R*), the concentration distribution of amendments in the matrix is more uniform at each depth resulting in a smaller residual rate, i.e., (R*-1) × 100%. If the injection wells have been installed unreasonably (e.g., a too large spacing), the continuous injection time is an effective controllable parameter to compensate for this defect. Moreover, a controlled slow-release system can maintain a more stable concentration distribution in the fracture than continuous injection and periodic injection systems, giving a longer residence time. Overall, this work is expected to provide some interesting guidelines for the design of multi-point injection in the fracturing low-permeability sites to enhance the remediation of contaminated soil.

Transient migration behavior of VOC vapor in layered unsaturated soils subjected to multiple time-dependent point pollution sources: Analytical study.

Predicting the migration behavior of volatile organic compounds (VOCs) vapor is essential for the remediation of subsurface contamination such as soil vapor extraction. Previous analytical prediction models of VOCs migration are mostly limited to constant-concentration nonpoint sources in homogeneous soil. Thus, this study presents a novel analytical model for two-dimensional transport of VOCs vapor subjected to multiple time-dependent point sources involving transient diffusion, sorption and degradation in layered unsaturated soils. Two representative time-dependent sources, i.e., an instantaneous source and a finite pulse source, are used to describe the short-term and long-term leakage. Results reveal that soil heterogeneity can cause pollution accumulation, especially in low-diffusivity capillary fringe. The assumption of an equivalent plane source from multiple point sources would significantly overestimate the vapor concentration and the contaminated range. The previous single point source model is no longer inapplicable when the relative distance and/or the release interval between sources is small, giving a strong interaction between multiple sources. Moreover, a faster vapor degradation rate or a higher groundwater level will reduce the area of vapor plume linearly. Hence, close attention should be paid to the time-variation characteristics of multiple sources, the vapor degradation and the groundwater level fluctuation in practice to facilitate soil remediation. The proposed model is a promising tool for addressing the above issue.

Design method of a modified layered aerobic waste landfill divided by coarse material.

To overcome the weaknesses of traditional landfills, a modified aerobic landfill concept with intermediate covers of coarse material between waste layers functioning as facilities of drainage and aeration has been proposed recently. In this study, a one-dimensional coupled model, including aerobic biodegradation, oxygen diffusion, and advection, is proposed to describe oxygen distribution in this new type of landfill. Homotopy analysis method and perturbation method are applied to solve this model at passive aeration and active aeration, respectively. The model has six input variables, that is, oxygen diffusion coefficient, gas permeability, maximum oxygen consumption rate, layer thickness of waste, and injection pressure and extraction pressure. A combination of their typical values gives rise to over 700,000 scenarios which can be calculated by the proposed solution. The coupled effect of the above variables on oxygen migration is quantitatively investigated, followed by an estimation formula of the minimum oxygen concentration in waste layer. The maximum waste layer thickness is defined as a function of other variables for a given aeration target of oxygen volume concentration larger than 5%. A generalized design method of waste layer thickness, injection pressure, and extraction pressure is then developed for the newly proposed modified layered aerobic landfill, which can promote its popularization and application.

Control and estimation of maximum gas pressure below landfill cover with horizontal gas wells: Analytical study.

Horizontal spacing of horizontal extraction gas wells can be designed to achieve a 90% pumping rate of the total generated landfill gas (LFG) from given waste properties (viz: gas permeability, landfill gas generation and non-homogeneity with depth), cover characteristics and vacuum pressure. However, cover characteristics and vacuum pressure are also important design parameters and different combinations of them result in different distributions of gas pressure in the waste, some of which would induce problematic air intrusion while others might pose threat to cover stability. This paper uses the maximum gas pressure directly below cover to distinguish these combinations, and provides the first study of the effects of the above parameters on potential outcomes. The ability of the overlying cover to resist LFG emission from the landfilled waste is suggested not to exceed a critical value, otherwise the maximum gas pressure below it would become at least 1 kPa larger than atmospheric pressure. A design formula for this critical value is proposed with respect to typical values of waste properties, vacuum pressure and the buried depth of horizontal wells in wide ranges. Together with consideration of recovery efficiency, the proposed method can be used to design a horizontal extraction gas collection system and a cover system in better working condition, and to evaluate the maximum gas pressure below cover. These applications are illustrated by a worked example.

Effect of LCRS clogging on leachate recirculation and landfill slope stability.

Vertical wells are commonly used for recirculating leachate into a landfill which can offer significant environmental and economic benefits. However, in some cases, the leachate collection and removal system (LCRS) at the bottom is overloaded and clogged due to biological and chemical processes. This results in a relatively high leachate level which could pose a threat to landfill slope stability. This study develops a three-dimensional landfill slope model with vertical recirculation wells and then investigates the effect of LCRS clogging on leachate recirculation and slope stability in terms of leachate saturation, pore water pressure, and factor of safety (FS) of a landfill slope. The results show that with an increase in clogging level that is characterized by an increased leachate level, the pore water pressure below the well injection screen is significantly increased by leachate recirculation, giving rise to a decreased slope FS value. In such conditions, the landfill slope formed by highly anisotropic waste is more likely to suffer instability. To prevent this kind of slope failure, a safe injection pressure of vertical recirculation wells is proposed for a wide range of parameter combinations involving waste anisotropy, clogging level, and the setback distance from the slope surface. This design guideline can be used to control the injection pressure in leachate recirculation applications and contributes to a better understanding of the slope stability of a bioreactor landfill.

Design of horizontal landfill gas collection wells in non-homogeneous landfills.

Considering exponential decreases in gas permeability and gas generation of waste with depth, a two-dimensional analytical model is developed to describe the landfill gas (LFG) recovery using horizontal wells. This model is used to simulate more than 680,000 scenarios involving typical values of waste properties, cover characteristics and design parameters for horizontal wells (seven variables in total). The coupled effect of these seven variables on air intrusion and the gas recovery efficiency of horizontal wells are investigated. It is shown that all the variables examined, except for the two variables defining waste non-homogeneity, could be integrated into three dimensionless variables. The horizontal spacing and buried depth of horizontal wells are examined as a function of cover characteristics, waste properties, and vacuum pressure to allow the development of a generalized design method for horizontal wells. An upper limit of horizontal well spacing is defined (for an 85% recovery rate) and a simple formula is provided which can be used to estimate the corresponding level of air intrusion. The upper limit spacing is shown to be affected by the non-homogeneity in gas permeability of waste, cover characteristics, and buried well depth. Using a worked example, the proposed method is shown to be capable of estimating air intrusion into existing horizontal gas collection wells and to optimize the design of horizontal wells considering waste non-homogeneity.

Numerical model of aerobic bioreactor landfill considering aerobic-anaerobic condition and bio-stable zone development.

Aeration by airflow technology is a reliable method to accelerate waste biodegradation and stabilization and hence shorten the aftercare period of a landfill. To simulate hydro-biochemical behaviors in this type of landfills, this study develops a model coupling multi-phase flow, multi-component transport and aerobic-anaerobic biodegradation using a computational fluid dynamics (CFD) method. The uniqueness of the model is that it can well describe the evolution of aerobic zone, anaerobic zone, and temperature during aeration and evaluate aeration efficiency considering aerobic and anaerobic biodegradation processes. After being verified using existing in situ and laboratory test results, the model is then employed to reveal the bio-stable zone development, aerobic biochemical reactions around vertical well (VW), and anaerobic reactions away from VW. With an increase in the initial organic matter content (0.1 to 0.4), the bio-stable zone expands at a decreasing speed but with all the horizontal ranges larger than 17 m after an intermittent aeration for 1000 days. When waste intrinsic permeability is equal or greater than 10-11 m2, aeration using a low pressure between 4 and 8 kPa is appropriate. The aeration efficiency would be underestimated if anaerobic biodegradation is neglected because products of anaerobic biodegradation would be oxidized more easily. A horizontal spacing of 17 m is suggested for aeration VWs with a vertical spacing of 10 m for screens. Since a lower aeration frequency can give greater aeration efficiency, a 20-day aeration/20-day leachate recirculation scenario is recommended considering the maximum temperature over a reasonable range. For wet landfills with low temperature, the proportion of aeration can be increased to 0.67 (20-day aeration/10-day leachate recirculation) or an even higher value.

Recovery response of vertical gas wells in non-homogeneous landfills.

A two-dimensional axisymmetric and normalized analytical model for landfill gas (LFG) migration around a vertical well is developed. The vertical gas permeability and LFG generation rate of waste are assumed to be subject to exponential decreases with depth. Using a general analytical solution, over 500,000 scenarios involving a combination of typical control variables (viz: cover properties, waste properties, vacuum pressure, well radius and spacing) are modelled. A quantitative analysis of the coupled effects of these control variables on LFG recovery rate indicates that the recovery response could be captured by: (a) three dimensionless variables (denoted as cover resistance, pump capacity, and well spacing parameters), and (b) two constants defining the decreases in gas permeability and LFG generation of waste with depth. For example, if the LFG generation rate of the waste at the top is doubled, a two times increase in the vacuum pressure with other parameters being equal would give a same gas recovery rate, as well as simultaneously doubling the thickness and gas permeability of the cover. The recovery efficiency of a vertical well with a low permeability cover is examined as a function of cover resistance and pump capacity, and design charts are presented that may be used to optimize gas recovery by adjusting cover properties and vacuum pressure. The proposed model makes it possible to consider the waste non-homogeneity in the design process, and the results contribute to a preliminary design of a cover and vertical LFG collection systems.

Design of vertical landfill gas collection wells considering non-homogeneity with depth.

Considering variable gas permeability and gas generation of waste with depth, different combinations of cover properties, vacuum pressure, and horizontal spacing of vertical wells giving rise to a 90% gas recovery rate are identified for typical waste properties. The effects of passive and active gas collection on horizontal well spacing are quantified. The normalized well spacing for 90% recovery is examined as a function of the cover resistance and the vacuum pump capacity. Design charts dependent on changes in gas permeability and gas generation rate with depth are presented to aid in the design of vertical gas wells. It is demonstrated that the non-homogeneity in gas permeability of waste is of great importance. For a conservative design, uncertainty in the non-homogeneity in gas permeability should be overestimated while uncertainty with respect to the non-homogeneity in LFG generation should generally be underestimated. The use of the proposed method for the design of the spacing of vertical gas wells in a situation with waste non-homogeneity is illustrated by a practical example.

CFD modeling of hydro-biochemical behavior of MSW subjected to leachate recirculation.

The most commonly used method of operating landfills more sustainably is to promote rapid biodegradation and stabilization of municipal solid waste (MSW) by leachate recirculation. The present study is an application of computational fluid dynamics (CFD) to the 3D modeling of leachate recirculation in bioreactor landfills using vertical wells. The objective is to model and investigate the hydrodynamic and biochemical behavior of MSW subject to leachate recirculation. The results indicate that the maximum recirculated leachate volume can be reached when vertical wells are set at the upper middle part of a landfill (H W/H T = 0.4), and increasing the screen length can be more helpful in enlarging the influence radius than increasing the well length (an increase in H S/H W from 0.4 to 0.6 results in an increase in influence radius from 6.5 to 7.7 m). The time to reach steady state of leachate recirculation decreases with the increase in pressure head; however, the time for leachate to drain away increases with the increase in pressure head. It also showed that methanogenic biomass inoculum of 1.0 kg/m3 can accelerate the volatile fatty acid depletion and increase the peak depletion rate to 2.7 × 10-6 kg/m3/s. The degradation-induced void change parameter exerts an influence on the processes of MSW biodegradation because a smaller parameter value results in a greater increase in void space.

A gas flow model for layered landfills with vertical extraction wells.

This paper developed a two-dimensional axisymmetric analytical model for layered landfills with vertical wells. The model uses a horizontal layered structure to describe the waste non-homogeneity with depth in gas generation, permeability and temperature. The governing equations in the cylindrical coordinate system were transformed to dimensionless forms and solved using a method of eigenfunction expansion. After verification, the effects of different well boundary conditions and gas extraction systems on recovery efficiency were investigated. A dimensionless double-layer system, consisting of a cover and a waste layer, was also explored. The results show that a constant vacuum pressure boundary condition can be enough to describe a perforated pipe surrounded by drainage gravel with a reasonable value of well radius, such as half the radius of gravel fill. Also, the 7 independent variables (one marked with an asterisk is dimensionless) of a double-layer system can be integrated into 3 dimensionless ones: Cover permeability Kv1∗/(Vertical gas permeability of waste Kv2∗×Cover thickness h1∗),-Vacuum pressure pw×PatmKv2∗/(µRgT2×Gas generation rate of waste s2) and ln(Well radius rw∗)/(Anisotropy degree of waste k2∗). The integration is based on the inherent mechanism of this flow system with certain simplification. The effects of these variables are then quantitatively characterized for a better understanding of gas recovery efficiency. Same recovery efficiency can be achieved with different variable combinations. For example, increasing h1∗ (such as doubling it) has the same effect with decreasing Kv1∗ (such as halving it). Along with the reduction of variables by half, the integration can facilitate the preliminary design, and is a small but important advance in the consideration of MSW non-homogeneity.

Unsaturated flow parameters of municipal solid waste.

Leachate pollution/recirculation and landfill gas emission are the major environmental concerns in municipal solid waste (MSW) landfills. A good understanding and prediction of MSW unsaturated properties are critical for the design of piping systems and the control of these problems within landfills. This paper reviews the recent studies of unsaturated properties of MSW, including experimental methods, theoretical models and corresponding model parameters. For experimental methods, the sample size is a common and significant limitation and large test apparatuses (e.g., >80cm in diameter) are generally required and valuable. The theoretical models for MSW also have some limitations due to the changes in waste composition and particle size distribution caused by biodegradation. Thus, the available data of intrinsic permeabilities, water retention curves, relative permeabilities and anisotropy of MSW were summarized to investigate the influences of porosity, waste composition and particle size distribution. A series of estimation methods were subsequently proposed to determine the parameters of water retention curve like θLm, θLr, nv and α. The other parameters such as the pore connectivity term (l) and the degree of anisotropy (k) were significantly lacking data, thus only their relationships with porosity were proposed. The results show that it is possible to define the second order effects caused by variations in porosity, waste composition and particle size distribution. However, the estimation methods still need more experimental data for improvement, especially their dependence on waste composition and particle size distribution.